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FIG. 1: (Color online) Panels a-b: The I — V characteristics of 
the A'^ = 10 channel IRLM calculated in the framework of the 
perturbative approach of Ref.Q], for inelastic rate Fin — and 
SFo, respectively. For finite Coulomb interaction U, I exhibits 
a power law asymptotics at large V. The corresponding expo- 
nent shown in panel c, essentially independent of Fi„. Panel 
d shows the mechanism leading to the decay of the current: 
spectral weight is transferred outside the voltage window. 


Comment on "Twofold Advance in the Theo- 
retical Understanding of Far-Prom-Equilibrium 
Properties of Interacting Nanostructures" 

Boulat, Saleur and Schmitteckert (BSS)[l[ report results 
on the full I—V characteristics of the interacting resonant 
level model (IRLM) exhibiting region with unexpected 
negative differential conductance (NDC). Using time- 
dependent density matrix renormalization group comple- 
mented with the exact solution performed at a special 
point (the self-dual point) in the parameter space BSS 
have shown that at nonzero Coulomb interaction U the 
current flowing through the impurity level (IL) exhibits 
a power-law asymptotics as a function of large applied 
bias voltage. Similar conclusion was earlier reached by 
Doyon[2]. Even though their results are solid and sup- 
ported by both analytic and numeric arguments, BSS 
concluded that "the NDC at large voltage seems a truly 
nonperturbative behavior, with unclear physical origin". 
On the contrary, the remarkable physics of NDC can be 
explained by simple physical arguments and, as we shall 
show in this Comment, in certain circumstances can be 
calculated in the framework of perturbation theory. 

As it was shown by Vladar and the present authors [1], 
by increasing the number of screening channels TV in the 


IRLM, the turning point of the hybridization exponent 
(cf. the self-dual point for iV = 2) can be pushed down 
to the perturbative regime. The scattering matrix el- 
ements were evaluated up to order U^, in the leading 
order of RG method. Using the scattering matrix ele- 
ments we can derive rate equations (RE) to determine 
the occupancy of the IL as well as the current / flowing 
through it. The RE approach takes only the sequential 
tunneling (ST) through the IL into account neglecting 
coherent co-tunneling (CT) processes, in contrast to the 
calculation of BSS. Nevertheless, for small hybridization 
F ST is dominant thus our results are reliable in that 
regime. To check that, we have repeated the calculation 
with the inclusion of a large inelastic rate Ti„ = SFq to 
de-phase electrons on the IL therefore suppressing CT 
(Fq being the bare value of F). In panels a-b of Fig. [2 
the I — V curves are shown for Fi„ = and = 5Fo, re- 
spectively. / clearly exhibits a power-law asymptotics 
at large V. As shown in panel c, the corresponding ex- 
ponent practically does not depend on the value of Fi„ , 
supporting the validity of ST approximation. It is also 
shown in panel c that the exponent of / coincide with the 
equilibrium exponent of F(a;) ~ Fo[a;/D]~^""'"^"' , where 
u is the dimensionless Coulomb interaction u = Ugo, 
go being the density of states (DoS) of conduction elec- 
trons per channel. Panel d illustrates the mechanism 
leading to NDC for u = 0.1: by plotting the IL DoS 
properly weighted with F-s (cf. the combination which 
enters the expression of the current) it is clear that spec- 
tral weight is transformed outside the voltage window. 
In the ST approximation the exponent of / coincide 
with the equilibrium exponent of F. For large enough 
V, the expression of / can be simplified dramatically, 
/ - T{eV) = TLieV) + TuieV). The saturation value 
of / is given by F which gets renormalized in presence 
of U but in non-equilibrium situation the flow is termi- 
nated at w ~ eV, thus I ~ r{eV) = ro[eV/D]-'^''+^"\ 
In general, our result does not rely on perturbative ar- 
guments: having the exact exponent of F from a reliable 
equilibrium calculation (e.g. numerical RG) one should 
be able to accurately describe the asymptotics of / in the 
ST regime. Note that in the calculation of BSS Fg was 
not very small therefore the CT was not negligible. That 
is why we found only a qualitative agreement between 
their exponents of / and those of F extracted from NRG. 

We conclude that in the ST regime of IRLM the NDC 
is a result of renormalization of F. This research was 
supported by Hungarian Grants OTKA No. T048782. 
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